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Abstract 

We comparetTtfie version 5 Microwave Limb Sounder (MLS) aboard the Upper Atmosphere 
Research Satellite (UARS), version 3 Polar Ozone and Aerosol Measurement-Ill (POAM-III) 
aboard the French satellite SPOT-IV, version 6.0 Stratospheric Aerosol and Gas Experiment II 
(SAGE-II) aboard the Earth Radiation Budget Satellite, and NASA ER-2 aircraft measurements 
made in the northern hemisphere in January-February 2000 during the SAGE III Ozone Loss 
and Validation Experiment (SOLVE). This study addresses one of the key scientific objectives of 
the SOLVE campaign, namely, to validate multi-platform satellite measurements made in the 
polar stratosphere during winter. This intercomparison was performed using a traditional 
correlative analysis (TCA) and a trajectory hunting technique (THT). Launching backward and 
forward trajectories from the points of measurement, the THT identifies air parcels sampled at 
least twice within a prescribed match criterion during the course of 5 days. We found that the 
ozone measurements made by these four instruments agree most of the time within ±10% in the 
stratosphere up to 1400 K (~35 km). The water vapor measurements from POAM-III and the 
ER-2 Harvard Lyman-a hygrometer and JPL laser hygrometer agree to within ±0.5 ppmv (or 
about ±10%) in the lower stratosphere above 380 K. The MLS and ER-2 CIO measurements 
agree within their error bars for the TCA. The MLS and ER-2 nitric acid measurements near 
17-20 km altitude agree within their uncertainties most of the time with a hint of a positive 
offset by MLS according to the TCA. We also applied the AER box model constrained by the 
ER-2 measurements for analysis of the CIO and HNO3 measurements using the THT. We found 
that: (1) the model values of CIO are smaller by about 0. 3-0.4 (0.2) ppbv below (above) 400 K 
than those by MLS and (2) the HN0 3 comparison shows a positive offset of MLS values by ~1 
and 1-2 ppbv below 400 K and near 450 K, respectively. It is hard to quantify the HN0 3 offset 
in the 400-440 K range because of the high sensitivity of nitric acid to the PSC schemes. Our 
study shows that, with some limitations (like HN0 3 comparison under PSC conditions), the 
THT is a more powerful tool for validation studies than the TCA, making conclusions of the 
comparison statistically more robust. 
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1- Introduction. 

The SAGE-III Ozone Loss and Validation Exper- 
iment (SOLVE) was a measurement campaign de- 1 
signed to study the processes controlling ozone behav- 
ior at high and mid-latitudes (see http://cloudl.arc. 
nasa.gov/soIve/ for details). SOLVE measurements 
were made mostly in the Arctic region during Novem- 
ber 1999 - March 2000 using the NASA DC-8 and ER- 
2 aircraft together with balloons and ground-based 
instruments. This mission was also intended to ac- 
quire correlative data necessary for validation of the 
SAGE-III measurements, which will be used to assess 
the behavior of global ozone and its possible trends. 
However, the launch of the SAGE-III instrument was 
postponed until 2001. Here, comparison of the ER-2 
measurements during the SOLVE campaign against 
available satellite instruments (POAM-III, MLS, and 
SAGE-II) is performed. This information can be used 
to benefit SAGE-EE validation later. 

Validation of any new instrument is a must before 
its products can be used for scientific studies. Tra- 
ditional correlative analysis (TCA) compares similar 
products obtained by new and other well-established 
platforms which are co-located in time and space as 
closely as possible. Such an approach is particularly 
attractive for comparison of the remote sensing and 
in situ data, since typically in situ measurements are 
more precise and accurate. To assess the quality of 
the satellite measurements, they are typically vali- 
dated with in situ measurements from ozonesondes, 
balloon, or aircraft (see the special issues of J . Geo- 
phys. Res., 94 (D6), 8335-8446, 1989; 101(06), 9539- 
10,476, 1996; 102(T> 19), 23,591-23,672, 1997 devoted 
to validation of SAGE-II, UARS, and POAM-II data, 
respectively). However, a limited altitude range and 
relatively small amount of matches (particularly for 
the occultation instruments with only 28-30 profiles 
per day [e.g., Lu et al , 1997a, b]) with in situ mea- 
surements could hamper a statistically significant val- 
idation of new satellite platforms. 

Recently, several new techniques using trajectories 
have been applied to improve validation of satellite 
data. Pierce et al. [1994] created “synoptic” maps 
of UARS Halogen Occultation Experiment (HALOE) 
measurements in order to statistically improve the 
comparison between HALOE sunrise and sunset mea- 
surements. Morris et al [1995] applied trajectory 
mapping to validate spatially-distant UARS MLS and 
HALOE data. The reverse-domain-filling technique 
[e.g., Sutton et al, 1994] was used to create uni- 
formly gridded satellite data by initializing trajecto- 
ries at a regular grid and then assigning them val- 
ues of the satellite measurements using backward tra- 
jectories and their encounters with satellite observa- 
tions. Bacmeister et al [1999] used the Lagrangian 
approach to map ER-2 flights on November 2 and 4, 


1994 and CRISTA data after November 4, 1994 to a 
GMT noon on November 5, 1994 for comparison. Lin - 
genfelser et al [1999] compared ER-2 and HALOE 
ozone measurements in the low^er stratosphere using 
trajectories. Morris et al. [2000] showed that tra- 
jectory mapping is a very effective tool for compar- 
ing one sparse data set with one dense satellite data 
set (HALOE and MLS O3), comparing two sparse 
data sets (HALOE and SAGE-II 0 3 ), and estimating 
instrument precision using MLS H 2 0 measurements. 
Von der Gathen et al [1995] and Rex et al [1998] 
applied the Match technique to find air parcels sam- 
pled twice by ozonesondes and to calculate ozone loss 
rates for the matched parcels. 

For validation of short-lived species, a box model 
is required in order to account for photochemical 
changes along matched trajectories. Pierce et al 
[1997] applied photochemical calculations along tra- 
jectories to compare ER-2 and HALOE measurements 
and model calculations of radical species against ER- 
2 measurements. Danilin et al [2000] used the AER 
box model along matched air trajectories in December 
1992 and obtained reasonable agreement between the 
calculated and measured behaviors of CIO, CIONO2, 
HNO3, and aerosol extinction at 780 cm -1 during this 
episode. 

The goal of this study is not to provide a thorough 
validation study for the MLS, POAM-III, and SAGE 
II instruments during their whole operation periods, 
but rather to perform a multi-platform data compar- 
ison between each other and with the ER-2 measure- 
ments for a particular period in January- February 
2000. In order to achieve this goal, we will use both 
the traditional correlative analysis (TCA) and trajec- 
tory hunting technique (THT). The structure of our 
paper is the following: Section 2 considers the main 
features of the instruments and the episodes consid- 
ered. Section 3 briefly describes the methods used. 
Sections 4, 5, 6, and 7 compare ozone, water vapor, 
chlorine monoxide, and nitric acid measurements, re- 
spectively. Finally, section 8 summarizes the main 
findings of our study. 

2. Period Considered and 
Characteristics of the Instruments 
Used. 

2.1. Period Considered. 

Figure 1 shows latitude coverage by the MLS 
(blue), POAM-III (red), SAGE-II (black), HALOE 
(yellow), and ER-2 (green) instruments from January 
10 to March 20, 2000. The ER-2 aircraft arrived in 
Kiruna (Sweden, 68°N, 21°N) on January 14, 2000. 
It made 6 flights during its first deployment (on Jan- 
uary 20, 23, 27, and 31 and February 2 and 3) and 5 
flights during its second deployment (on February 26 
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and March 5, 7, 11, and 12). It returned to the USA 
on March 16. To support SOLVE, POAM-III changed 
its routine schedule when it looks every other day ei- 
ther in the northern (sunrise) or southern (sunset) 
hemisphere. The POAM-III provided sunrise mea- 
surements in the northern hemisphere almost every 
day (except February 14, 2000), with 14 profiles daily 
most of the time. Additionally, the MLS obtained 
data from February 2 to 13, 2000 and from March 27 
to 29, 2000 (the last period is not shown in Figure 1). 
To maximize the probability of obtaining useful mea- 
surements in the sunlit polar vortex, MLS scans were 
performed only poleward of 50°N on the “day” side 
of the orbit, providing several tens of profiles of O3, 
HN0 3 , and CIO daily [Santee et al, 2000b]. 

Two other occultation instruments (SAGE-II and 
HALOE) operated on their routine basis. SAGE- 
II provided sunrise measurements poleward of 50° N 
from January 29 to February 8, 2000. HALOE 
looked poleward of 50°N from February 12 to 27 (sun- 
rise) and from March 7 to 24 (sunset). Since the 
HALOE latitude overlapped barely with that of MLS 
on February 12-13 and briefly crossed the path of the 
ER-2 flights on March 11 and 16, very few matches are 
anticipated for the MLS/HALOE and ER-2/HALOE 
pairs. The poor statistics of the comparison for these 
pairs precludes us from using the HALOE data in our 
analysis. We focus our study on the available mea- 
surements during the January 23 - February 18, 2000 
period, which is characterized by the relatively dense 
coverage by the satellite and aircraft measurements. 

2.2. Instruments Used. 

Table 1 summarizes the instruments and their prin- 
cipal parameters used in this study. 

ER-2. The NASA ER-2 aircraft has 17 instru- 
ments aboard, of which the following five are used 
in this study for comparison with satellite measure- 
ments: dual-beam UV-absorption ozone photometer 
[Proffitt et al., 1989], the Harvard Lyman-a hygrome- 
ter [Weinstock et al , 1994], the JPL laser spectrome- 
ter [May, 1998], the Harvard N02-C10-C10N02-Br0 
instrument [Stimpfle et al , 1999], and the Caltech 
CIMS instrument [McKinney et al. , 2001]. Here we 
briefly describe the main parameters and principles of 
measurements of these instruments. 

Ozone is measured by the ozone photometer, which 
uses a light from_ 254 nm lamp. This light passes 
through two identical sample chambers, one with am- 
bient air and the other with the same air with com- 
pletely removed ozone. The absorption cross section 
of ozone is well-known and large at 254 nm. Thus, 
the difference between the detected signals from two 
chambers allows accurate determination of the ozone 
concentration in the ambient air. At a measurement 
frequency of 1 Hz, the minimum detection limit of 


ozone is about 1.5xl0 10 molecules/cm 3 (or 8 ppbv at 
20 km) with a total uncertainty of several percent. 

The Harvard Lyman-a photofragment fluorescence 
hygrometer measures water vapor aboard the ER-2 
aircraft. Ambient flow is ram fed through the nose 
of the aircraft. To optimize accuracy, the core of the 
flow is picked off and throttled rapidly through the 
water vapor detection axis (for a full description of the 
instrument, see Weinstock et al. [1994]). Reported 
accuracy of the measurements is ±5% with a potential 
offset of ±0.1 ppmv [Hintsa et al , 1999]. The l-a 
precision of the measurements is typically ±0.1 ppmv 
for a 10-second integration time. 

Water vapor is also measured by the JPL laser hy- 
grometer [May, 1998]. This is a single-channel, near- 
infrared tunable diode laser spectrometer operating 
near 1.37 micrometers wavelength with a multipass 
optical cell in the Herriott configuration. To ensure no 
contamination, the open-path optical cell is mounted 
external to the boundary layer of the right wing su- 
perpod of the ER-2. This instrument is calibrated in 
the laboratory with a known water vapor standard 
over the range of pressures experienced during flight, 
and compared with a chilled mirror frost-point hy- 
grometer. For SOLVE data, the one-sigma precision 
is typically better than ±0.05 ppmv over a 1.3-second 
integration time. The reported accuracy is ±5% at 
pressures <100 hPa, and ±8% at pressures between 
100 and 200 hPa [May, 1998]. 

Chlorine monoxide is measured by the Harvard 
N02-C10-C10N02-Br0 instrument, which comprises 
two separate instruments: a thermal dissociation, res- 
onance fluorescence instrument for detection of halo- 
gen radical and reservoir species [Stimpfle et al, 1999] 
and a laser-induced fluorescence instrument for the 
detection of NO2 [Perkins et al, 2001]. Within the 
halogen system, CIO is measured by reaction with in- 
jected NO to form Cl atoms, followed by resonance 
fluorescence detection of Cl atoms at 118.9 nm. The 
instrument sensitivity is calibrated in the laboratory 
with known Cl atom densities and normalized to the 
Rayleigh scattering signal. The CIO measurements 
are acquired with a 35 sec temporal resolution with 
an uncertainty (lcr) and detection limit of ±17% and 
3 pptv, respectively. 

HNO3 is measured by chemical ionization mass 
spectrometry (CIMS). The instrument inlet selec- 
tively samples either gas or particles in the ambient 
air using a modified virtual impactor technique. Typ- 
ically, aerosol and gas phases are sampled alternately 
during flight for periods of 3 minutes each. The sam- 
ple then passes through a flow tube at ~290 K to 
evaporate condensed HNO3 from aerosols. HNO3 is 
ionized by chemical reaction with the precursor ion, 
CF 3 0": CF3O- + HN0 3 HF-NO^ + CF 2 0. Ions 
are then directed through an aperture into the vac- 
uum system, where, following mass selection, they are 
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detected by a channel electron multiplier. A constant, 
known amount of isotopically labeled HN0 3 is added 
to the sample, allowing continuous, simultaneous cal- 
ibration and measurement. For the data used here 
(gas phase measurements from the January-February 
ER-2 deployment) the precision is about ±0.75 ppbv 
(lcr) for a 7-s integration period. The accuracy is the 
greater of ±25% and ±1 ppbv. 

MLS, The microwave limb sounding technique and 
the MLS instrument are described in detail by Wa- 
ters [1993] and Barath et al. [1993], respectively. 
The MLS instrument began acquiring millimeter- 
wavelength emission measurements in late September 
1991. MLS measurements were made daily or near- 
daily until 1995, after which data coverage became 
progressively sparser; in July 1999 the MLS instru- 
ment was placed in standby mode and was not pow- 
ered up again until February-March 2000, when the 
205-GHz radiometer was operated to provide mea- 
surements of CIO, 0 3j and HN0 3 in the northern 
hemisphere in conjunction with the SOLVE cam- 
paign. Santee et al. [2000b] present these measure- 
ments and describe the ways in which the operational 
strategy employed when these data were collected dif- 
fered from the previous mode. The year 2000 data 
may have small shifts (typically of order a few per- 
cent) in comparison to normal v.5 operational re- 
trievals (with the 63 GHz radiometer providing tan- 
gent pressure and temperature information), but the 
2000 data are still considered useful for the type of 
analyses performed here. In this study we analyze 
the CIO, 0 3 , and IINO3 measurements obtained dur- 
ing the February 2000 observing period. 

A special issue of Journal Geophysical Research 
( 101(D6 ), 9539-10,476, 1996) was devoted to the 
UARS validation and discussed the MLS version 3 
(v.3) results in detaiL Subsequently, version 4 data 
were released; additional information on version 4 
ozone data was provided in WMO [1998]. Version 5 
(v.5) MLS data have recently become available. In 
v.5, quantities are retrieved on every UARS surface 
(six surfaces per decade in pressure, as opposed to 
three in previous MLS data sets), although the true 
vertical resolution of the data has not doubled. Over- 
all, the data quality and the vertical range of relia- 
bility has been improved in v.5, especially for ozone 
in the lower stratosphere. Estimated vertical reso- 
lutions and total uncertainties of the individual mea- 
surements of ozone, nitric acid, and chlorine monoxide 
are about 4 km, 6 km, and 4 km and about 0.4 ppmv, 

2 ppbv, and 0.4 ppbv, respectively, at the altitudes of 
interest in this study. A detailed description of the 
MLS v.5 data processing algorithm and validation of 
the various v.5 data products will be given elsewhere 
[ N . Livesey et al, manuscript in preparation]. Infor- 
mation about the quality of the v.5 data is also avail- 
able from the MLS web site (http://mls.jpl.nasa.gov). 


POAM-III. The POAM-III instrument is a nine 
channel visible/near-infrared photometer which mea- 
sures O3, H2O, NO2, and aerosol extinction using 
the solar occultation technique. It was launched on 
March 23, 1998 into a Sun-synchronous 98.7° incli- 
nation orbit with a period of 101.5 min and is still 
in operation. Occultation measurements provide 14 
sunrise and 14 sunset profiles daily with a vertical 
resolution of about 1-2 km (depending on species and 
altitude) and a longitudinal distance of 25.7° between 
successive measurements. Because of the orbit incli- 
nation, the measurements cover the 55°N-71°N and 
62°-88 0 S latitude bands providing sunrise and sunset 
profiles in the northern and southern hemispheres, re- 
spectively. The center wavelengths of the POAM-III 
channels are located at 354, 439.6, 442.2, 603, 761.3, 
779, 922.4, 935.9, and 1018 nm. A detailed descrip- 
tion of the POAM-III and its early validation results 
are presented in Lucke et al [1999] and at the POAM 
web site (http://opt.nrl.navy.mil/POAM). Also, the 
special POAM-II validation section of JGR [102, 
no.D19, 23,591-23,672, 1997] could be useful for the 
POAM-III data analysis, since these two instruments 
are similar. 

SAGE-II. The SAGE-II instrument was launched 
into a 57° inclination orbit aboard the Earth Radia- 
tion Budget Satellite on October 5, 1984 [Russell and 
McCormick , 1989] and is still operational. SAGE- 
II employs the solar occultation technique and mea- 
sures O3, NO2, H2O, and aerosol extinctions with 
a vertical resolution of better than 1 km. This in- 
strument measures attenuated solar radiation at 7 
wavelengths centered near 0.385, 0.448, 0.453, 0.525, 
0.60, 0.94, and 1.02 microns. Details of the SAGE- 
II retrieval algorithm and validation of the SAGE-II 
measurements are given elsewhere [e.g., Chu et al , 
1989; JGR f 94 , D6, 8335-8446, 1989; http://www- 
sage2.larc.nasa.gov]. We use the latest version 6.0 of 
the SAGE-II data in this study. 

While papers detailing the v.6.0 algorithm and val- 
idation of the data products are in preparation, an 
overview of the changes implemented in version 6.0 
can be found at http://www-sage2.larc.nasa.gov/cd- 
rom. The primary motivation for the version 6.0 de- 
velopment was to understand and correct the long 
standing bias in the ozone profiles between SAGE-II 
and ozonesondes in the 15 to 20 km region [WMO, 
1998]. Version 6.0 development was also focused on 
improving the overall quality of the SAGE II ozone 
and aerosol products. The primary improvements 
were to the transmission algorithm, in the way of im- 
proved altitude registration, and to the methodology 
of calculating and propagating the known sources of 
error. Compared to the previous version, v.5. 96, this 
new version of ozone has finer vertical resolution, con- 
siderably smaller error estimates, reduced aerosol ar- 
tifacts, and a greater vertical range extending into the 
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mid troposphere. 

3. Analysis Techniques. 

3.1. Traditional Correlative Analysis. 

TCA is a conventional method used for validation 
of atmospheric measurements during the last several 
decades. TCA finds nearly-coincident profiles mea- 
sured by the different platforms that satisfy a pre- 
scribed match criterion. The averaged profiles mea- 
sured by different instruments are compared and their 
differences (in absolute units or in %) are analyzed in 
order to find averaged offsets between different in- 
struments. This approach is used, for example, to 
compare UARS and SAGE-II measurements with ra- 
diosonde or ground-based measurements \J. Geophys . 
Res., 94 (D6), 8335-8446, 1989; 101{ D6), 9539-10,476, 
1996; WMO, 1998], The main drawback of the TCA 
is the small number of matches. In order to improve 
the statistics of the comparison, more relaxed match 
criteria can be used. However, care is needed in order 
to avoid compromising the comparison due to using 
an overly-relaxed criterion. For example, ozone be- 
comes a relatively short-lived gas above 35 km and 
its diurnal variability can spoil a comparison of the 
nadir versus occultation measurements if a temporal 
match criterion of several hours is used. 

Additionally, comparison of zonally averaged data 
can be made. In this case, meteorological variabil- 
ity and small scale effects are mostly removed and 
the statistical significance of such a comparison is im- 
proved. However, the very large temporal and spa- 
tial difference between the measurements used in this 
approach is a shortcoming. For comparison of short- 
lived species (like NO2 or N2O5), the solar zenith an- 
gle of measurements is important in order to avoid 
diurnal modulation of the species of interest. Despite 
all these caveats and limitations, the TCA remains 
the most trusted and widely used method of multi- 
platform data validation and special validation cam- 
paigns are usually planned for a new platform (like 
EOS Aura). 

3.2. Trajectory Hunting Technique. 

THT identifies air parcels sampled at least twice 
by the same or different platforms and compares mea- 
surements along the matched trajectories [Danilin et 
al, 2000]. There are four stages in applying the THT 
for validation studies. At the first stage, backward 
and forward trajectories are calculated from the loca- 
tions of measurements of interest. In this study, only 
5-day backward and forward trajectories are used for 
all pairs. We calculated 670 backward and forward 
trajectories that originated from the ER-2^ points at 
p<110 hPa (i.e. above ~15 km) for the ER-2 flights 
on January 23, 27, 31, and February 2 and 3, 2000. 


Since the frequency of measurements varies for the 
different instruments aboard ER-2, for the sake of 
convenience we use the ER-2 merged files provided 
by R.J. Salawitch with 10 s averaging. However, the 
trajectories were calculated from the ER-2 points dis- 
tanced with an interval of 2 min along the flight track. 
Comparing SAGE-II and MLS ozone data, we cal- 
culated 4044 backward and forward trajectories be- 
tween 15 and 35 km with a 1 km step originated from 
the SAGE-II profiles made during the February 1-15, 
2000 period. For comparison of POAM-III and MLS 
ozone measurements, 5130 backward and forward tra- 
jectories are calculated between 9 and 35 km altitude 
with a 1 km vertical step originated from the POAM- 
III profiles measured during the February 1-18 pe- 
riod. We did not calculate trajectories above 35 km 
(~1400 K), since ozone has a short photochemical life- 
time (less than 1 day) above this level and could not 
be considered as a passive tracer. Also, mixing with 
ambient air becomes rapid in the mesosphere and up- 
per stratosphere [Shepherd et al, 2000], thus limiting 
the effectiveness of the THT there. Since MLS does 
not provide reliable measurements for p>100 hPa, we 
also did not consider trajectories below this level. 

The diabatic trajectories used in this study were 
computed using 3-D winds derived from the temper- 
ature and geopotential height fields from the U.K. 
Meteorological Office assimilation scheme [Swinbank 
and O'Neill, 1994]. The mean meridional wind and 
the vertical wind (mean and eddy components) were 
calculated using the thermodynamic and continuity 
equations in the same manner as in Smith and Ly- 
jak [1985]. The net diabatic heating rates used in the 
thermodynamic equation were calculated as described 
in Gille and Lyjak [1986]. The three-dimensional wind 
at the locations and times required by the trajectory 
calculation was obtained by linearly interpolating the 
daily wind fields in space and time. Our selective 
comparison shows that our trajectories are very close 
to those provided by the Goddard Automailer. 

At the second stage of the THT, we check whether 
each trajectory launched from the locations of the first 
platform measurements passes within a prescribed 
temporal-spatial distance from the other platform 
measurements. If yes, the hunting for this trajectory 
is successful and it is a subject for further analysis; if 
not, we drop this trajectory from our study. It is bet- 
ter to launch trajectories from the relatively sparse 
measurements (ER-2 or POAM-III) and to hunt for 
the more frequent measurements (like MLS) rather 
than vice versa, since this approach requires less com- 
puter time and memory while providing the same re- 
sults. 

The third stage is devoted to the interpolation in 
the vertical coordinate of the matched measurements. 
This interpolation is required only for the targeted 
measurements (i.e., MLS), since the values at the ini- 
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tial points of the trajectories are known* For the re- 
sults shown below we used a linear interpolation in 
the log-pressure scale. However, our sensitivity analy- 
sis shows that a choice of the vertical coordinate (e.g., 
pressure vs log-pressure or 9 vs pressure) does not no- 
ticeably affect the results. Once the matched points 
are known, the same interpolation code is applied for 
all species sampled by both platforms. 

At the final, fourth stage, grouping and statisti- 
cal analysis of the matched measurements are per- 
formed. We define “grouping” as a procedure that 
bins all matched data as a function of any vertical 
coordinate (potential temperature, pressure or alti- 
tude). For example, below we bin all satellite/satellite 
matched measurements with a step of 50 K and 100 K 
from 350 K to 1000 K and above 1000 K, respectively. 
We increased the grouping step above 1000 K in or- 
der to get comparable statistics with the data below 
1000 K, since a vertical step in km per 100 K in po- 
tential temperature is smaller in the middle than in 
the lower stratosphere. For the ER-2/MLS and ER- 
2/POAM-III pairs we apply a constant vertical step 
of 20 K. 

This brief description of the THT shows that TCA 
matches are only a subset of the THT matches, when 
a match is obtained near initial points of trajecto- 
ries. Thus, the THT is a statistically more powerful 
tool than the TCA for validating atmospheric mea- 
surements. THT is also a more cost-efficient way to 
carry out validation campaigns, since it allows for ob- 
taining as much useful information as possible from 
independent measurements in the background atmo- 
sphere in addition to specially deployed platforms. In 
our paper we show comparisons for the same pair of 
instruments using both the TCA and THT in order 
to increase confidence in the results shown. 

For all but the MLS/SAGE-II results shown below, 
we use the match criterion of (Atime<2h, Alatitude<2 £> , 
Alongitude<2°). For the Kiruna latitude, the spatial 
difference of 2° in latitude and longitude is translated 
into 237 km. For the MLS/SAGE-II ozone measure- 
ments, we apply the match criterion of (Atime<8h, 
Alatitude<2°, Alongitude<3°), since for a shorter 
Atime no matches are found for the period consid- 
ered using the TCA. Below, for the sake of simplicity, 
we write the match criterion as, for example, (2 h, 2° ) 
2°) omitting the words “Atime<”, “Alatitude<”, and 
“Alongitude<”. In order to facilitate analysis of the 
results presented, we will show in some cases the dif- 
ference between two instrument measurements both 
in absolute (ppmv or ppbv) and relative (%) units. 

4. Comparison of the Ozone 
Measurements. 

The figures comparing MLS with satellite (POAM- 
III and SAGE-II) and ER-2 measurements are shown 


with a vertical resolution of ~ 1-1.5 km or ~0.5 km, 
respectively, while the vertical resolution of the v.5 
MLS data is about 4-6 km. This caveat should be 
kept in mind. We decided to show a finer vertical 
resolution linearly interpolating the MLS data in the 
log(pressure) coordinate to the data by POAM-III, 
SAGE-II, and ER-2. The results shown below look 
quite similar and conclusions of our study are not af- 
fected when we use a vertical step of ~4 km after 
degrading the POAM-III and SAGE-II vertical reso- 
lution to the MLS vertical grid. For the ER-2/MLS 
pair the MLS vertical resolution provides only one 
level for comparison. 

4.1. MLS vs POAM-III. 

We found only two MLS/P OAM-III profiles satis- 
fying the chosen match criterion of (2 h, 2°, 2°), which 
are shown in Table 2. We averaged the matched MLS 
and POAM-III profiles separately and then show their 
difference by black lines in Figure 2. Launching 5-day 
trajectories from the POAM-ni points, we found a to- 
tal of 3113 matches satisfying the same match crite- 
rion. The red line in Figure 2 shows the difference be- 
tween the POAM-III and MLS ozone measurements, 
when all these matches are grouped with the 50 K and 
100 K steps below and above 1000 K, respectively. 

Both methods show that POAM-III ozone values 
are smaller than those of MLS almost everywhere be- 
low 1400 K. This difference ranges from 0 to -0.6 ppmv 
(or from 0 to 12%) for the THT and from 0 to -1 
ppmv (or from 0 to 19%) for the TCA. Results for the 
TCA could be suspicious because of the small statis- 
tics (i.e., only 2 individual profiles). However, the 
THT results have about 100-300 individual matches 
per level shown and should not suffer from possible 
large discrepancies in the individual ozone measure- 
ments. The effect of photochemical changes along the 
matched trajectories is unlikely to be an issue for this 
comparison, since the number of matches for forward 
(1605 matches) and backward (1446 matches) trajec- 
tories is about the same (see Table 3 for statistics 
of the THT). This factor indeed should be consid- 
ered if we are in the region of intensive photochemi- 
cal ozone destruction. In this case, POAM-III ozone 
values will be reduced (increased) by the photochem- 
istry for the forward (backward) trajectories origi- 
nating from the points of POAM-III measurements. 
However, detailed analysis of the photochemical ozone 
loss rate shows a modest ozone depletion in the begin- 
ning of February due to the lack of sunlight [Hoppel 
et ai, 2001]. Also, our box model calculations for the 
ER-2/MLS ozone comparison show only a small pho- 
tochemical ozone loss in the lower stratosphere (see 
sections 4.3 and 6 for details). All these arguments 
convince us that the average difference between the 
POAM-III and MLS ozone measurements is real. 
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4.2. MLS vs SAGE-II. 

Table 2 shows that because of the geometry of 
the SAGE-II and MLS measurements, there are no 
matches within a 6 h interval. To get at least a min- 
imum statistic for the MLS/SAGE-II comparison us- 
ing the TCA, we relaxed the match criterion from (2h, 
2 a , 2°) to (8h, 2°, 3°). For this criterion, we found 
5 matched profiles. Again, we averaged them sepa- 
rately for each instrument and depict the difference 
between their averaged values by black lines in Fig- 
ure 3. The difference obtained by the THT is shown 
by red lines and the number of the matches found is 
depicted near the right vertical axis. We found 3397 
matches for the 4044 trajectories launched. 

The obvious feature of this Figure is the good 
agreement between the v.5 MLS and v.6.0 SAGE- 

II measurements. The difference never exceeded the 
range of il0%(±5%), being within ±5% (±3%) at 
most levels for the TCA(THT). The largest differ- 
ence of 7-8% is obtained near 1000 K (~28 km). It is 
likely that this difference arises from the small number 
of statistics obtained for the TCA, when the effects 
of possible large differences for individual profiles is 
weakly mitigated for the averaged profiles. Indeed, 
we saw a large difference (MLS - SAGE-II > 1 ppmv) 
at 26-32 km for the first of five matched ozone profiles 
listed in Table 2, which was propagated in Figure 3. 
The statistically significant difference in the range 
of [-2%:-4%] and [-5%:-9%] is obtained in the mid- 
dle stratosphere (900 K to 1200 K) for the THT and 
TCA, respectively. This difference is smaller than 5% 
of the MLS positive offset for the v.3 MLS and v.5. 93 
SAGE-II comparison above 30 hPa reported by Cun- 
noli et al. [1996] and Froidevaux et ai [1996], thus 
showing an improvement obtained in the latest ver- 
sions of these data sets. The very good agreement be- 
tween MLS and SAGE-II throughout the stratosphere 
shows no evidence of the large biases seen in earlier 
versions of the SAG El- II algorithms. Our comparisons 
of the ozone measurements for the MLS/SAGEMI and 
MLS/POAM-III pairs are consistent with the results 
of Manney et al [2001] showing that the agreement 
between v.5 MLS and six other instruments (includ- 
ing v.6.0 SAGE-II and v.5 POAM-II) is better than 
0.25 ppmv in the stratosphere. Summarizing, such a 
good agreement between MLS, SAGE-II, and POAM- 

III validates the use of the latest versions of MLS, 
POAM-III and SAGE-II ozone measurements for sci- 
entific purposes. 

4.3. MLS vs ER-2. 

There are two ER-2/MLS matches satisfying the 
match criterion of (2h, 2°, 2°) obtained during the 
February 2, 2000 flight using the TCA. During this 
flight, ER-2 flew at its cruise altitude between the 
two consecutive MLS profiles taken at llh37m and 


llh38m UT (see also Table 2). The ER-2 stayed 
14 min within the chosen vicinity from each of the 
coincident MLS profiles. Figure 4 depicts the ER-2 
and MLS ozone values obtained during these matches. 
The MLS ozone values are provided at 100, 68, and 
46 hPa, while ER-2 sampled at 60-65 hPa. The ER-2 
and MLS ozone measurements agree within their un- 
certainty for these two episodes. However, a possible 
positive offset in the MLS data of 0. 2-0.3 ppmv is sus- 
pected. When the match criterion was relaxed to (3 h, 
2°, 5°), six matches were found for the February 2, 
2000 ER-2 flight (not shown). For all these matches, 
the ER-2 and MLS ozone measurements agree within 
their uncertainties with even a closer agreement than 
shown in Figure 4. In such comparisons, one needs 
to keep in mind the huge difference in the individual 
volume sampled (400x400x4 km 3 for MLS and much 
less than 1 km 3 for ER-2, which does not sample much 
of the altitude region sensed by MLS). This and the 
small number of matches makes it difficult to assess 
any systematic difference (outside the error bars) be- 
tween the two sets of measurements, especially for the 
TCA method. 

The THT allows better quantification of the dif- 
ference between the MLS and ER-2 ozone measure- 
ments, since more matches are found covering a wider 
vertical range. We found 525 matches for the 387 
trajectories launched from the ER-2 points during its 
flights on January 31, February 2 and 3, 2000. Other 
ER-2 flights are irrelevant for comparison with MLS, 
since they occurred more than 5 days before the start 
or after the end of the MLS measurements in Febru- 
ary. Figure 5 shows that ER-2 ozone measurements 
are smaller than those from MLS by 0.2-0. 3 ppmv 
in the 360 K - 460 K range. These results are sta- 
tistically more robust than those of Figure 4. Part 
of the possible small discrepancy between MLS and 
ER-2 data may be linked to small offsets obtained 
for the single-radiometer mode retrievals used after 
mid- 1997, in comparison with the multi- radiometer 
retrievals used prior to this with the 63 GHz data 
available for tangent pressure and temperature re- 
trievals. Indeed, tests performed by the MLS team 
indicate that the single-radiometer retrievals tend to 
overestimate the results from the full-blown retrievals 
by 0.1 to 0.2 ppmv at high northern latitudes. Our 
box model calculations initialized with the ER-2 mea- 
surements showed a modest averaged photochemical 
ozone depletion of 0.005 ppmv during the averaged 
temporal distance of 1.8 days between the matched 
ER-2 and MLS measurements (see Section 6 for de- 
tails). Thus, the possible photochemical changes 
along the matched air parcels did not affect the re- 
sults of the ER-2/MLS ozone comparison and can be 
ignored. 


4.4. POAM-III vs ER-2. 

Using the TCA, we found two ER-2/P0AM-III 
matches on February 2 and January 23, when ER-2 
was cruising and descending to Kiruna, respectively. 
These matches are shown in Figure 6, which depict 
10 second averaged ER-2 data and part of the corre- 
sponding coincident POAM-III profile. The ER-2 was 
10 min and 52 min within (2h, 2°, 2°) from the coinci- 
dent POAM-III profile. While the ER-2 and POAM- 
III ozone measurements show good agreement within 
their uncertainties, one can suspect a possible small 
positive bias in the POAM-III measurements. In or- 
der to check this assumption, we relaxed the match 
criterion to (3h, 2°, 5°) and found 7 matches (not 
shown). However, they showed a similar behavior 
of the POAM-III and ER-2 measurements shown in 
Figure 6. Applying the THT to the POAM-III/ER- 
2 pair, we found 405 matches, which are shown in 
Figure 7 after grouping them with a vertical step of 
20 K. The agreement between ER-2 and POAM-III 
is good (better than 8%), especially taking into ac- 
count a huge difference in individual volume sampling 
of the POAM-III (200x60x1 km 3 ) and ER-2 ( much 
less than 1 km 3 ). 

Lumpe et al [2001] also compared ozone measure- 
ments made by ER-2 and POAM-III during SOLVE. 
In general, their results are very similar to ours, 
confirming good agreement between the ER-2 and 
POAM-III ozone measurements. For example, they 
showed an agreement between ER-2 and POAM-III 
ozone measurements within ±10% in the 350 K - 
470 K range using both the vortex-averaged and tra- 
jectory matching techniques. 

5. Comparison of the ER-2 and 
POAM-III H 2 0 Measurements. 

POAM-III measures H 2 0 using the 935.9 and 922.4 
nm channels. Previous comparisons of v.1.4 POAM- 
III and v.19 HALOE measurements showed a high 
bias of about 15% in the stratosphere in the POAM- 
III measurements [Lucke et al., 1999]. Recently, 
Bevilacqua et al. [2001] compared v.3 POAM-III and 
v.19 HALOE H 2 0 measurements and found that the 
POAM-III measurements are higher by <10% in the 
20-40 km range in the northern hemisphere. The ER- 
2 water vapor measurements have been made during 
many aircraft campaigns and involved in the recent 
SPARC H 2 0 assessment [SPARC, 2001]. 

Figure 8 shows the ER-2/POAM-III comparison 
for water vapor. Unlike the ozone case, there are two 
instruments aboard ER-2 which measure H 2 0 (see 
section 2 for their details). During the February 2, 
2000 flight, both ER-2 instruments showed very com- 
pact water vapor values of about 5 ppmv at 62 hPa, 
almost overlapping each other. POAM-III measured 
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slightly higher water content, which could be consid- 
ered consistent with the ER-2 measurements within 
uncertainties of the measurements. During the de- 
scent on January 23, 2000, both the Harvard and JPL 
instruments showed an increase of H 2 0 which is also 
hinted in the POAM-III measurements. However, 
individual water vapor measurements by POAM-III 
can be noisy in the lower stratosphere/upper tropo- 
sphere, especially if sunspots and clouds are present 
(see Bevilacqua et al [2001] and Nedoluha et al 
[2000] for details). For example, the POAM-III H 2 0 
profile shown in Figure 8b is strongly affected by a 
PSC presence at this location according to the sharply 
increased aerosol extinction detected by the POAM- 
III aerosol channels. 

Water vapor behaves as an inert tracer in the 
stratosphere, unless temperature drops below the ice 
frost point causing removal of H 2 0 from the gas 
phase. However, based on the analysis of tempera- 
ture along the trajectories used, such cold tempera- 
tures did not occur during the period considered here. 
Thus, neglect of possible microphysical changes is eas- 
ily justified for the use of THT for H 2 0. We found 402 
and 354 matches for Harvard and JPL hygrometers, 
respectively. The number of matches is smaller for the 
JPL instrument because of some gaps in its data dur- 
ing three ER-2 flights (000131, 000202, and 000203). 
Figure 9 shows the difference between the H 2 0 mea- 
surements by the ER-2 instruments and POAM-III, 
which changes from ±8% (-3%) at 370 K to -8% (- 
12%) at 470 K for the JPL (Harvard) hygrometer. 
Most of the time at cruise altitude, the difference be- 
tween the Harvard and JPL hygrometers lies within 
the reported uncertainties [Hintsa et al, 1999]. How- 
ever, there is a pressure-dependent systematic bias 
between the two instruments that is greatest at 100 
to 200 hPa pressure, which has been attributed to 
the JPL laser hygrometer [SPARC, 2001]. As shown 
in Figure 9, this bias is about 0.4 ppmv at 390 K re- 
ducing to ~0.2 ppmv at 490 K (JPL values are always 
larger) . 

The SOLVE POAM-III/ER-2 H 2 0 comparison in 
our paper shows some inconsistency with the previous 
POAM-III/HALOE and ER-2/HALOE comparisons 
[SPARC, 2001]. For example, our study shows that 
v.3 POAM-III water vapor measurements are higher 
by 5-10% and 10-15% than the Harvard hygrometer 
and v.19 HALOE data, respectively, at 15-20 km al- 
titude. This implies that the Harvard data should 
be higher by 5-10% than the v.19 HALOE measure- 
ments. However, the SPARC [2001] study shows that 
the Harvard hygrometer data are 20% higher than 
v.19 HALOE measurements. This inconsistency is 
discussed elsewhere [Bevilacqua et al, 2001] and its 
further analysis is outside the scope of this paper. 
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6* Comparison of the ER-2 and MLS 
CIO Measurements. 

6.1. Analysis using the TCA. 

Previously version 3 MLS CIO measurements were 
validated against balloon, aircraft and ground-based 
measurements [Waters et al, 1996]. It was found that 
the v.3 MLS CIO measurements agree with available 
correlative measurements within their combined un- 
certainties. Subsequently, v.4 MLS CIO data were 
compared with the 204 GHz measurements from the 
millimeter-wave atmospheric sounder (MAS) on the 
spac^ shuttle [Feisi et al. , 2000]. The agreement was 
well within the combined error bars over a pressure 
range of 0.4-40 hPa. In v.5 MLS, use of a finer ver- 
tical grid in the retrievals allows better definition of 
the peak and results in generally smoother profiles 
(although the true vertical resolution of the v.5 CIO 
measurements, 4-5 km at the altitudes shown here, 
is coarser than the retrieval grid). A positive bias of 
0.1 ppbv in the lower stratosphere is known to exist 
in the v.5 CIO data (based on averages of nighttime 
data over the first full year of the mission) [N. Livesey 
et 'al.y manuscript in preparation]. Here we compare 
the v.5 MLS data with SOLVE measurements. 

Figure 10 compares CIO measurements made by 
the Harvard N02-C10-C10N0 2 -Br0 instrument aboard 
the ER-2 (green dots) and MLS (blue lines). For the 
two ER-2/MLS matches shown, the measurements 
by these instruments agree within their uncertainties. 
During these matches, chlorine was activated, with 
values between 1 and 2 ppbv in the 70-45 hPa range. 
Such a small number of matches can be understood 
from Figure 1, which shows that only two ER-2 flights 
(on February 2 and 3) are available for direct com- 
parison with MLS measurements. Unfortunately, the 
February 3 flight was very short and remote from the 
MLS points by at least 10 h, making comparison of 
the measurements of short-lived CIO impossible for 
this flight using the TCA. To increase the number 
of coincident ER-2 and MLS CIO measurements, we 
relaxed the match criterion to (2h, 2°, 5°). In this 
c ase, we found 6 matches (not shown here). In all 
these matches ER-2 and MLS agreed within their er- 
ror bars. However, ER-2 values of CIO were smaller 
by at least 0.1 ppbv, consistent with the known high 
bias in the MLS CIO data at these altitudes. 

The agreement between remote sensing values of 
CIO in a box 400x400x4 km 3 and very accurate in 
situ ER-2 measurements is good. It is clear that lower 
values of MLS CIO in v.5 compared with their values 
from the previous versions mitigate a concern about 
their positive bias raised previously by modelers [e.g., 
Chipperfield et al, 1996; Lutman et al, 1997] and 
mentioned by Waters et al. [1996]. Since chlorine 
monoxide is a key ozone depleting species, accurate 
measurements of CIO and their agreement with mod- 


els are crucial for our understanding of the changes in 
the global ozone layer. 

6.2. Analysis using the THT. 

Because of the short photochemical lifetime of CIO 
in the stratosphere, the THT requires a photochem- 
ical box model for the matched parcels. We use the 
AER box model [Danilin et al , 2000] constrained by 
the ER-2 measurements for analysis of the CIO mea- 
surements. The same model runs are used for analysis 
of the HNO 3 measurements presented in section 7.2. 
For the forward trajectories originated from the ER- 
2 points, our model is initialized using the ER-2 O 3 , 
NO, N0 2 , HN0 3 , NO y , CIO, C1 2 0 2 , HC1, CIONO2, 
H 2 0 (JPL hygrometer values), CH 4 , and N 2 0 values. 
Initial concentrations of other chlorine and nitrogen 
species, which are not measured by the ER-2 (namely 
HOC1, Cl 2 , OCIO, HONO, N 2 0 5 , and HN0 4 ), are 
determined from the ER-2 measurements of NOy and 
Cl y using the partitioning among these species ac- 
cording to the AER 2-D model at 67°N in February 
[ Weisenstein et al, 1998]. However, the initial con- 
centrations of these six species are small, contribut- 
ing only several percent to the NO y and Cl y amounts 
and are not important for our comparison of the CIO 
and HNO 3 measurements. In order to justify the last 
statement, we performed a sensitivity model run with 
zero initial concentrations of these six species and ob- 
tained the same results for CIO and HNO 3 as shown 
below. Some details of nitrogen species partitioning 
between NOy and HNO 3 in gas and condensed phases 
do not affect our model CIO results and are given in 
section 7.2. Initial total inorganic bromine Br y is de- 
termined using the N 2 0-Br y correlation based on the 
analysis of Wamsley et al [1998] and is about 18 
pptv at the ER-2 cruise altitude. The initial parti- 
tioning among bromine species is determined accord- 
ing to the AER 2-D model. The initial aerosol sur- 
face area density (SAD) is taken according to the Fo- 
cused Cavity Aerosol Spectrometer (FCAS) [Jonsson 
et al, 1995] measurements. We discuss later why the 
choice of the SAD measured by FCAS or Multiple- 
Angle Aerosol Spectrometer Probe (MASP) [Baum- 
gardner et al, 1996] has little impact on our results. 
Our model has an option to use either nitric acid tri- 
hydrate (NAT) [Hanson and Mauersberger, 1988] or 
supercooled ternary solution (STS) [Tabazadeh et al, 
1994] PSC schemes. Details of our model treatment 
of PSCs are given elsewhere [Danilin et al, 2000]. 
Briefly, in the NAT scheme, the PSC SAD appears 
at T<T NA t in addition to the sulfate aerosol SAD 
(which was constant in this scheme) and is propor- 
tional to the amount of HNO3 condensed. In the STS 
scheme, the aersol SAD increases depending on the 
amount of HNO3 and H 2 0 condensed at low temper- 
ature. Also, the reaction probablities are different on 
the NAT and STS surfaces and are taken according 
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to Sander et al., 2000]. One also should keep in mind 
that the amount of HNO3 condensed in the NAT and 
STS scheme is quite different, especially in vicinity 
of Tyv^r- Below we applied both schemes to analyze 
CIO and HNO3 measurements. 

It was suggested that the UKMO assimilated tem- 
perature has a warm bias in the polar lower strato- 
sphere during winter [e.g., Manney et al, 1996]. In- 
deed, we found that the UKMO temperature was 
higher on average by 1.2 K compared with the ER-2 
measurements in the initial points of the 450 forward 
trajectories. To offset this warm bias, the forward 
trajectory temperatures were forced to match ER-2 
temperatures in the initial points. One should keep 
in mind that despite this procedure, the assimilated 
temperature in other points of trajectories could still 
deviate from its “real” values. This procedure is stan- 
dard in model analysis of ER-2 measurements during 
polar winters [e.g., Kawa et al t 1997], allowing for 
considerable reduction of uncertainties of this crucial 
parameter (±0.3 K for ER-2 vs several K for UKMO). 

We performed model calculations for the 450 for- 
ward trajectories originated from the locations of ER- 

2 measurements. We did not make model runs for the 
backward trajectories from the ER-2 locations, since 
only MLS O3, CIO, and HNO3 values are available 
in the initial MLS points (i.e., where model calcu- 
lations start), thus introducing additional uncertain- 
ties of the model initialization due to unknown ini- 
tial partitioning of the chlorine and nitrogen species. 
Also, we notice some inconsistency between the MLS 
HNO3 and ER-2 NO y measurements, since the MLS 
HNO3 values often (in 14 of 21 cases) were larger 
than the ER-2 NO y for forward matches, thus mak- 
ing model initialization difficult. Due to the timing 
of the ER-2 and MLS measurements (see Figure 1), 
the forward trajectories provide the dominant part of 
all ER-2/MLS matches (450 of 525, or ~86%). In- 
deed, since MLS started its measurements on Febru- 
ary 2, only about 0.5 and 2 days of the MLS data are 
match-available for the backward trajectories origi- 
nated from the ER-2 locations during the February 2 
and 3 flights, respectively. On the other hand, 5 and 

3 days of the MLS measurements are match- available 
for the forward trajectories during these two flights 
and the January 31 flight, respectively. 

The pink line in Figure 11 shows the vertical pro- 
file of the difference between the MLS-measured and 
model-calculated CIO values using the NAT PSC 
scheme. This difference has maximum values of about 
-0.4 ppbv near 390 K and smaller values of about - 
0.2 ppbv above 400 K. For the STS scheme, the re- 
sults (black line) are almost the same as for the NAT 
scheme. The main reason for this is the very low ini- 
tial values of C10N0 2 (<0.05 ppbv) and HOC1 (a few 
tens of pptv). Thus, despite the difference in the re- 
action probabilities and SADs for the NAT and STS 


schemes, the C10N0 2 + HC1 -> Cl 2 + HN0 3 and 
HOC1 + HC1 -» CI2 + H 2 0 heterogeneous reactions 
are shut down, precluding any additional formation 
of active chlorine. 

It is important to notice that the error bars in Fig- 
ure 11 show the standard error of the mean differ- 
ences (as in all previous Figures) and do not account 
for systematic error of the MLS CIO measurements. 
If this error is also taken into account, the error bars 
in Figure 11 should be increased three- fold below 400 
K and two-fold above this level. Model initialization 
and calculations also introduce additional uncertain- 
ties. However, thorough analysis of this very impor- 
tant issue is outside the scope of this paper, requires 
numerous sensitivity calculations (as in Considine et 
al [1999] or Cohen et al. [2000]) and deserves a sep- 
arate study. 

7. Comparison of the ER-2 and MLS 
HNO3 Measurements. 

7.1. Analysis using the TCA. 

Designed primarily to measure stratospheric abun- 
dances of CIO, O3, and H 2 0, UARS MLS also mea- 
sures HNO3. Previous v.4 MLS measurements of 
HNO3 are discussed in detail by Santee et al [1999]. 
Recently, Santee et al [2000a] presented a first 
preliminary validation of the v.5 MLS HNO3 data 
with Atmospheric Trace Molecular Spectroscopy (AT- 
MOS) and UARS Cryogenic Limb Array Etalon Spec- 
trometer (CLAES) HNO3 measurements. The v.5 
MLS HNO3 values agree well with both ATMOS 
(v.3.1) and CLAES (v.9) data at most altitudes in 
the tropics. At mid-latitudes, the agreement between 
these platforms is generally good, except for the 46-22 
hPa range, where the MLS data are higher by up to 
35% and 50% compared with ATMOS and CLAES, 
respectively. Under conditions when HNO3 is en- 
hanced inside the winter polar vortices in regions of 
limited PSC activity, MLS values can exceed those 
of the infrared measurements by 15-60%. Detailed 
validation of the v.5 MLS HNO3 data is currently in 
progress [Santee et al, manuscript in preparation]. 
The CIMS is a new instrument aboard ER-2 and op- 
erated for the first time during the SOLVE campaign. 
CIMS validation studies are still under way [ McKin- 
ney et al, 2001]. 

Figure 12 shows the comparison of the matched 
ER-2 and MLS measurements on February 2, 2000. 
Because the CIMS instrument measures gaseous HNO3 
during approximately half of the flight time (during 
the other half it measures condensed-phase HNO3), 
and because of the longer integration time for the 
HNO3 data, fewer measurements of nitric acid are 
reported for a given ER-2 flight compared with mea- 
surements for ozone. Consequently, fewer ER-2 points 
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are shown in Figure 12 than in Figure 4. In order to 
increase the number of matches, we relaxed the match 
criterion to (3 h, 2°, 10 a ). If the match criterion of 
(2 h, 2°, 2°) is used, only the two matches in the top 
row of Figure 12 remain. It is difficult to quantify the 
difference based on the TCA results shown in Fig- 
ure 12 because of the poor statistics available and a 
narrow vertical range of 52-63 hPa covered by ER-2 
during its February 2 flight. Usually, the MLS and 
CIMS ER-2 HNO3 measurements agree within their 
uncertainties. On average, perhaps, one can say that 
the ER-2 values tend to be smaller than the MLS data 
(especially in panels (a) and (e)). However, panel (f) 
and in lesser extent panel (b) show the cases when 
ER-2 values are larger than those by MLS. 

7.2. Analysis using the THT. 

To improve the statistical significance of the MLS/ER- 
2 nitric acid comparison, we applied the THT with a 
match criterion of (2 h, 2°, 2°). We found a total of 
148 matches (127 and 21 matches for the forward and 
backward trajectories, respectively), less than for the 
ozone and CIO cases because of fewer available HNO3 
measurements by ER-2. The results of this analysis 
are shown by the red, blue and green lines in Fig- 
ure 13a for all, forward, and backward trajectories, 
respectively. 

The results shown in Figure 13a assume that gaseous 
nitric acid is a passive tracer along the matched tra- 
jectories. For a good passive tracer, the difference be- 
tween the forward and backward trajectories should 
be small, indicating no changes in the tracer concen- 
tration with time. However, the large difference (up 
to 4 ppbv at 430 K) between results for the forward 
and backward trajectories in Figure 13a clearly sig- 
nals that nitric acid is not a passive tracer above 
400 K. We investigated whether the changes in the 
HNO3 concentrations are caused by the following 
reasons: (a) photochemistry, (b) heterogeneous re- 
actions, (c) irreversible denitrification of the HNO3- 
containing PSC particles, and (d) temporary removal 
of nitric acid from gas to condensed phases and back 
via condensation to and evaporation from the PSC 
particles, respectively. 

(a) The photochemical changes of nitric acid are 
small in our study (less than few tens of pptv accord- 
ing to our model calculations below). However, this 
statement could be easily confirmed without any de- 
tailed model calculations using the fact that, on aver- 
age, the matched air parcels were illuminated about 
6 h. Assuming a HNO3 photolysis rate of the or- 
der of 10 -7 - 10“ 6 s -1 , one gets that about 0.2-2% 

(or ~0.01-0.1 ppbv) of HNO3 can be destroyed. This 
value is consistent with our model calculations and 
much smaller than the ER-2/MLS differences shown 
in Figures 12-13. 


(b) The effects of heterogeneous reactions on sul- 
fate aerosol and PSCs are also small for HNO3. The 
main reason for this is very small initial values of ni- 
trogen (NO, N0 2 , N 2 0 5 and C10N0 2 ) species during 
this period, which could be eventually converted to 
HNO3 via heterogeneous reactions. All these species 
have concentrations smaller than several tens of pptv, 
which could produce at most 0.1 ppbv of additional 
HNO3. The last value is an order of magnitude 
smaller than the differences we saw in Figure 13a. 
We performed our model runs using the FCAS aerosol 
SAD. While usually the MASP SAD values are larger 
than those by FCAS, we do not anticipate that our 
comparisons of the HNO3 ER-2 /MLS measurements 
will be affected because of the small values of NO, 
N0 2 , N 2 0 5j and C10N0 2 . 

(c) In this study we assume that there is no 
irreversible denitrification caused by sedimentation 
of HN03-containing particles between the matched 
points. This assumption is supported by the fol- 
lowing facts. First, temperature along the matched 
air trajectories was always above the ice frost point, 
thus precluding formation large ice particles. Second, 
Tabazadeh et al [2000] showed that irreversible den- 
itrification could happen in air parcels that stayed 
for ~1 week at temperatures below T nat an d above 
Ti ce . In our study, the averaged temporal differ- 
ence between ER-2 and MLS points is 1.9 days, too 
short to denitrify the air even assuming that all this 
time temperature was below T nat- However, the 
recent discovery of the large 10-20 micron particles 
(so called “rocks”) [Fahey et a/., 2001], which could 
come from higher layers and reduce (increase) the 
HNO3 amount because of sedimentation (evapora- 
tion), complicate the justification of the above as- 
sumption about the lack of irreversible denitrification 
between the matched points. Indeed, some of these 
“rocks” were obtained by the forward NO y instrument 
during some periods of the January 31 and February 3 
flights [Fahey et al, 2001]. We hope that they would 
not change the HNO3 amount in the matched parcels 
noticeably because of their large sedimentation speed 
and short stay in our parcels. However, we do not as- 
sess possible effects of the “rocks” on our ER-2/MLS 
HNO3 comparison and acknowledge that their possi- 
ble impact could introduce additional uncertainties. 

(d) Nitric acid may experience rapid changes caused 
by condensation to and evaporation from PSC par- 
ticles in the polar lower stratosphere during winter. 
Qualitatively, the difference for the forward and back- 
ward trajectories in Figure 13a could be understood if 
one takes into account temperatures in the ER-2 and 
MLS matched points, which is illustrated for the for- 
ward trajectories in Figure 13c. The ER-2 flew under 
temperatures below the NAT threshold most of the 
time during the three flights considered. For such 
conditions, HN0 3 may be removed from the gas to 
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solid or liquid phase. On the other hand, on average, 
temperatures in the MLS points were warmer than 
those in the ER-2 points. As a result, less HN0 3 may 
be condensed onto PSCs in the MLS locations. Thus, 
for the forward trajectories originated from the ER-2 
points, one should anticipate a release of the HN0 3 
condensed into the gas phase as the parcels move to- 
wards the MLS points, thus making the apparent val- 
ues of the ER-2 HNO3 higher. For the backward tra- 
jectories originated from the ER-2 points (not shown 
in Figure 13c), the situation is opposite, making the 
apparent ER-2 HNO3 values smaller. Based on this 
analysis, one also can say that the large difference be- 
tween the THT HNO3 results for forward and back- 
ward trajectories indicates possible PSC events cap- 
tured during the period considered. 

To account for the microphysical and, less im- 
portant, photochemical changes of HNO3 along the 
matched forward trajectories, we use the AER box 
model with NAT or STS PSC schemes [Danilin et 
ai, 2000]. Some information of the model initializa- 
tion relevant for analysis of the CIO measurements 
was given in section 6.2. Here, we provide some de- 
tails of the model initialization that are important for 
comparison of the HNO3 measurements. Ideally, the 
initial concentrations of all nitrogen species should be 
derived from the ER-2 measurements. However, some 
nitrogen species (like N2O5 or HNO4) are not mea- 
sured or are not available simultaneously (like HNO3 
in gas and condensed phases). We use the following 
approach in order to initialize our model runs. If the 
difference between the NOj, measured by the back- 
ward inlet of the NO y instrument [Fahey et al, 2001] 
and the sum of all other nitrogen species is positive, 
we assigned the balance to nitric acid condensed in 
aerosol (HNO3). We also checked that the non-zero 
HNO3 values are consistent with low temperature 
in the ER-2 locations. Indeed, Ter 2 was less than 
T nat in most (^90%) such cases. For the remain- 
ing (~10%) cases, we kept the non-zero HNOf value 
despite relatively warm temperature in these points. 
Also, the ER-2 measurements of nitrogen species were 
not always internally consistent and the CIMS HN0 3 
values were larger than the ER-2 NO y data in 15 of 
127 matches. The reason of this inconsistency is dis- 
cussed elsewhere [McKinney et al, 2001]. In such 
cases, we initialized our model using the ER-2 mea- 
surements of HNO3 , N0 2j NO and CIONO2, assum- 
ing concentrations of all other nitrogen species equal 
to zero. Even in such cases, we avoid any scaling 
of initial HNO3 that may hamper our further com- 
parison. These details are important for the analysis 
below. 

Figure 13b shows that the THT analysis of the 
MLS and ER-2 HNO3 measurements is a very chal- 
lenging task even for the forward trajectories, when 
we have the best currently available constraints by 


the ER-2 data on the model initialization and a rel- 
atively large amount of matches compared with the 
backward trajectories case. The blue line with error 
bars is reproduced from Figure 13a, while the black 
and pink lines show the model results when the pho- 
tochemical and PSC-related changes in HNO3 are^ac- 
counted for the matched parcels using the STS and 
NAT schemes, respectively. The obvious feature of 
this figure is a large difference between the model re- 
sults with different PSC schemes. 

Again, qualitatively this difference could be under- 
stood in terms of the averaged temperature in initial 
and final points of the matched trajectories. Below 
400 K level, the results for the STS and NAT calcula- 
tions are the same, since the averaged temperature in 
the final points was above the NAT threshold. Above 
400 K level, the averaged temperature in the final 
points drops below T nat, thus causing a big differ- 
ence in the amount of HN0 3 condensed onto PSC 
particles in NAT and STS schemes. It is well-known 
that more HNO3 condenses in the NAT than STS 
scheme in the temperature range of [T;v at-^'^nat\ 
[Tabazadeh et al. , 1994; Carslaw et al, 1994]. This 
fact is also widely used in order to discriminate be- 
tween different PSC schemes using satellite HNO3 
measurements [e.g., Santee et al, 1998]. In our study, 
we see the maximum difference of ~3 ppbv between 
the NAT and STS schemes at 430 K level. 

The yellow line in Figure 13b shows the results of 
the model runs with the STS PSC scheme and zero 
initial HNOf. We perform this calculation in order to 
study the sensitivity of the ER-2/MLS HN 0 3 com- 
parison to this parameter, which was not measured 
by ER-2 for the matches shown (otherwise, we can 
not compare measurements of gas-phase HNO3). The 
caveat for this scenario is the neglect of the ER-2 NO y 
measurements for the model initialization (since the 
sum of all nitrogen species concentrations is different 
from the ER-2 NO y values), thus making this model 
run less trustworthy. It is obvious that the yellow 
line is very close to the blue line and far from the 
black line. Thus, the difference between results for 
the model STS run (black line) and the tracer sce- 
nario is caused by the assumed initial HNOf, which 
averaged values were equal to 0.66, 1.03, 0.28, 0.78, 
and 0.28 ppbv at 370, 390, 410, 430, and 450 K, re- 
spectively, for our model runs. The barely noticeable 
difference of 0.04 ppbv between the yellow and blue 
lines at 370 and 390 K is caused by photochemistry, 
since the effects of HNO3 condensation onto PSC par^ 
tides are removed. Indeed, the initial HNO3 content 
in aerosol is equal to zero by the assumption for this 
scenario, while the final HNOf values are also equal 
to zero because of the high temperature in the MLS 
points (see Figure 13c). 

Results shown in Figure 13 show that the THT 
suffers for the ER-2/MLS HNO3 measurement com- 
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parison in the polar lower stratosphere during win- 
ter because of the rapid HNO3 condensation to or 
evaporation from PSC particles. Obtaining an accu- 
rate quantitative estimate of the difference between 
the MLS and ER-2 measurements using the THT 
is further complicated by its high sensitivity to the 
PSC scheme used. Based on the presented results, 
we can conclude that the MLS HNO3 measurements 
are larger by 1-2 ppbv at 450 K and by ~1 ppbv at 
370-390 K than those from the ER-2. It is difficult 
to quantify the difference in the 400-440 K range be- 
cause of the high sensitivity of the HNO3 results to 
the PSC schemes (NAT vs STS). 

A big difference in the vertical and horizontal res- 
olution between MLS and ER-2 measurements is an- 
other important caveat, which should be kept in mind. 
These factors limit the usefulness of the THT in pro- 
viding robust conclusions for the ER-2/MLS HNO3 
comparison in the polar lower stratosphere during 
winter. However, we anticipate that the THT analy- 
sis should be more successful for the HNO3 measure- 
ments in any other regions of the stratosphere that 
are free from PSC activity. 

8. Summary and Conclusions. 

The SOLVE campaign provided a unique oppor- 
tunity to validate remote sensing satellite measure- 
ments (POAM-III, MLS, and SAGE-II) with the in 
situ ER-2 measurements in the polar winter strato- 
sphere. We compared these platform measurements 
during the January-February 2000 period. Also, we 
estimated only the random error when comparing dif- 
ferent pairs of instruments and did not account for 
systematic errors and possible biases. We also de- 
cided that uncertainties introduced by the trajectory 
calculations and model initialization and calculations 
are outside the scope of this paper and deserve a sepa- 
rate study. Analyzing the latest versions of the MLS 
(v.5), POAM-III (v.3), SAGE-II (v.6.0), and ER-2 
measurements, we conclude that: 

(1) the ozone measurements by MLS and SAGE- 
II are in excellent agreement (better than 5% and 
mostly within 2%), with MLS values slightly larger 
than those from SAGE-II. The POAM-III ozone mea- 
surements are up to 12% smaller than those from 
MLS. The MLS ozone values are larger by 6%-ll% 
than those from the ER-2 in the 400-480 K vertical 
range, but huge sampling volume differences make it 
difficult to provide significant confidence in offset de- 
terminations between MLS and ER-2 measurements. 
The POAM-III and ER-2 ozone measurements show 
better than 8% agreement in the 380-500 K vertical 
range with some evidence of a very small (about 5%) 
POAM-III positive offset (also mentioned in [Lumpe 
et a/., 2001]). These differences between the vari- 
ous techniques typically fall within the expected com- 


bined accuracies of the different data sets and other 
assumptions of the methods used, which we estimate 
to be of the order of 10%. Such small differences be- 
tween the ozone measurements made by these plat- 
forms are very encouraging for scientific applications 
of these ozone data sets; 

(2) the POAM-III water vapor measurements are 
in agreement with both the Harvard Lyman-a and 
JPL laser hygrometers within ±0.5 ppmv (or about 
±10%) with a hint of a higher offset as large as several 
tenths of 1 ppmv in the POAM-III data in the 360- 
480 K range; 

(3) according to the TCA, the MLS and ER-2 CIO 
measurements agree within their error bars and MLS 
shows a positive offset consistent with the 0.1 ppbv 
bias known to be present in MLS v.5 CIO data in the 
lower stratosphere. Model calculations constrained 
by the ER-2 measurements show that the MLS mea- 
surements are higher by 0.2 (0.4) ppbv above (below) 
430 K. However, the uncertainties introduced by the 
model calculations are not evaluated in this study and 
can be large. Also, the sampling volumes for the MLS 
and ER-2 measurements are very different; 

(4) the MLS and CIMS ER-2 HNO3 measurements 
are consistent within their uncertainties most of the 
time with some hint of a positive offset by MLS ac- 
cording to the TCA. PSC processing and the high 
sensitivity of the HNO3 to the choice of PSC scheme 
complicate quantification of their difference in the 
400-440 K range using the THT. However, at 450 K 
and below 400 K the positive offset of MLS is 1- 
2 ppbv and about 1 ppbv, respectively, according to 
the THT. The statistical significance of these values 
is not easy to assess with high confidence, however, 
given the model uncertainties related to the PSC mi- 
crophysics discussed in section 7 and the different vol- 
ume sampling by MLS and ER-2. 

Table 3 summarizes some important statistical 
data characterizing the THT for the period studied. 
For example, we found more than 3000 matches for 
the MLS/POAM-III and MLS/S AGE-II pairs using 
the 5-day trajectories compared with 2 and 5 profiles 
(or 50 and 125 matches assuming 25 matches per pro- 
file), respectively, satisfying the same match criteria 
as in the THT. Table 3 confirms the fact that more 
matches are obtained for more frequent measurements 
and more relaxed match criteria. 

Our study shows that the trajectory hunting tech- 
nique is an effective tool in validation of multi-platform 
measurements that provides more statistically robust 
conclusions than the TCA, due to the much larger 
number of matches obtained. Compared with the tra- 
jectory mapping technique [Morris et al, 2000], we 
believe that the THT presented here is more computa- 
tionally efficient, since no mapping of the dense mea- 
surements (like MLS) is performed. Also, the THT 
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results do not depend on a chosen time for which the 
mapping is performed. However, it is possible that 
these techniques provide similar results regarding off- 
sets of one instrument against another one. 

In the future, we would like to investigate how suc- 
cessful the THT could be in conjunction with a photo- 
chemical box model for analysis of short-lived species 
and to estimate model uncertainties for such compar- 
isons. The methodology presented here could be used 
for validation studies of future space-borne platforms 
(like SAGE-III or EOS Aura). 

Appendix: Error Analysis. 

The standard deviation of the differences between 
any two platform measurements is defined according 

to (1): 


SD = 


N 


N 


) 2 /(^-l) 


( 1 ) 


here N is the number of matches, A* is the difference 
for the i-th pair (e.g., Ai=0^ LS -0f- Oj4M “ //7 ) and 
A*/N is the mean difference between the two 
instrument measurements. The values of A; could 
be expressed in absolute (ppmv or ppbv) or relative 
(%) units. The standard deviation characterizes the 
spread of the distribution near the mean value (i.e., 
A) and is a measure of the combined random error of 
both instruments. Thus, it does not account for the 
systematic errors. The error bars shown in all figures 
presenting the THT results are determined according 
to (2): 

ERR = SD/VN , (2) 

showing that the error bars become smaller for larger 
N. The values of ERR represent the standard error of 
the mean differences and correspond to the 67% con- 
fidence level. The values of ERR should be doubled 
for the 95% confidence. 
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Table 1. Instruments Analyzed in This Study and Their Characteristics (Measured Parameters with Their 
Uncertainty ±1 a. Vertical Resolution Az, and Solar Zenith Angle of Measurements) 



Instrument 

Parameters Measured 

Az, km 

SZA 


MLS 

Oa(±0.3 ppmv), HN03(±l-2 ppbv), C10(±0.3-0.6 ppbv) 

4-6 d 

any 


POAM-III 

O3(±0.5 ppmv), H20(±lppmv), NO2(±0.5 ppbv) 

1 

90° 


SAGE-II 

0 3 (±1%), NO 2 (±10%), H 2 O(±0.5 ppmv) 

0.5 

90° 

— 

ER-2° 

0 3 (±0.05 ppmv),H20(±5% 6 ,±5% c ),HN0 3 (±0.75ppbv),C10(±17%) 

<0.1 

any 


a Partial list of the parameters measured by ER-2 is shown; POAM-III and SAGE-II also measures aerosol extinctions 
b Harvard Lyman-c* hygrometer; C JPL laser hygrometer 
d 4 km for O3 ajid CIO and 6 km for HNO3 
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Table 2. Time and Locations of the Matches Between the Instruments Shown According to the TCA Satisfying 
the Match Criterion of (2h, 2°, 2°) 


Pair 

Date 

Time (UT) and L( 

MLS/POAM-III 

February 2 

llh37m, 66.36°N, 

MLS/POAM-III 

February 3 

22h52m, 65.34°N, 

MLS/SAGE- II a 

February 7 

19h50m, 52.38°N, 

MLS/SAGE-IP 

February 7 

21h44m, 52.33°N, 

MLS/SAGE-II 

February 8 

10hl5m, 51.92°N, 

MLS/SAGE-IP 

February 8 

13h28m, 51.82°N, 

MLS/SAGE-IP 

February 9 

10hl7m, 51.14°N, 

MLS/ER-2 6 

February 2 

llh37m, 66.36°N, 

MLS/ER-2 6 

February 2 

llh38m, 63.08°N, 

POAM-III/ER-2 fr 

January 23 

13h37m, 65.34°N, 

POAM-III/ER-2 ( ' 

February 2 

12h03m, 66.31°N, 


50.55°E / 12h03m, 66.3PN, 48.74°E •; 

123.99°W / 23h33m, 66.44°N, 122.74°W 
85.19°W / 13h23m, 50.76°N, 87.91°W 
109.52°W / 14h59m, 50.71°N, 112.04°W 
55.88°E / 3h50m, 50.24°N, 54.90°E 
7.22°E / 7h02m, 50.12°N, 6.63°E 
51.00°E / 3h53m, 49.16°N, 52.81°E 
50.55°E / llhOOm, 64.50°N, 48.74°E 
54.59°E / llhl4m, 64.72°N, 54.65°E 
17.97°E / 12h03m, 65.38°N, 19.00°E 
48.74°E / 10h55m, 64.41°N, 49.31°E 


a The match criterion of (8h, 2°, 3°) is used for the MLS/SAGE-II pair 
b Mean values of ER-2 time and location during the match are shown 
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Table 3. Statistical Data for the THT with the Match Criterion of (2 h, 2°, 2°) Shows the Total Number of 
Relevant Trajectories N tra and Matches NJ 71 , Number of Matches for the Backward (Nj*) and Forward (NJ 1 ) 
Trajectories, Number of Matches per Trajectory Launched (N^/Nt™), and Mean Temporal Distance between the 
Matches (At, days). 


Pair 

Species 

N^ ra 

wr 

Nr 


Nr/N (ra 

At, d 

MLS/POAM-III 

0 3 

5130 

3051 

1446 

1605 

0.59 

2.4 

MLS/SAGE-IP 

0 3 

4044 

3397 

907 

2490 

0.84 

2.2 

MLS/ER-2 

0 3 

387 

525 

75 

450 

1.36 

1.8 

MLS/ER-2 

CIO 

387 

525 

75 

450 

1.36 

1.8 

MLS/ER-2 

hno 3 

387 

127 

21 

105 

0.33 

1.9 

ER-2/POAM-III 

0 3 

670 

405 

148 

257 

0.60 

2.3 

ER-2/POAM-III 

h 2 o 6 

670 

402 

147 

255 

0.60 

2.3 

ER-2/POAM-III 

H 2 O c 

670 

354 

132 

222 

0.53 

2.3 


a Match criterion of (8 h, 2°, 3°) was applied for this pair 
b Harvard Lyman-a hygrometer; C JPL laser hygrometer 




ER-2-green; MLS-blue; POAM3-red; HALO E-yellow; SAGE-ll-black 


JANUARY 


FEBRUARY 

2000 


MARCH 


Figure 1: Latitudinal coverage of the MLS (blue), POAM-III (red), ER-2 (green), SAGE-II 
(black), and HALOE (yellow) measurements during the January 10 - March 20, 2000 period. 
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Figure 2: Vertical profiles of the difference between POAM-III and MLS ozone measurements 
in ppmv (panel a) and in % (panel b) from the TCA (black curves) and THT (red curves). 
Numbers of matches for the THT are shown near the right vertical axis. Match criterion of 
(2 h, 2°, 2°) is used. The error bars shown represent only the random (not systematic) error 
determined in Appendix and correspond to ±1 a (or 67% confidence level). 
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Figure 3: The same as in Figure 2, but for the SAGE-II and MLS ozone measurements and a 
match criterion of (8 h, 2°, 3 n ). 
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Figure 4: Two matches between ER-2 and MLS ozone measurements for the ER-2 flight on 
February 2, 2000 according to TCA with the match criterion of (2 h, 2°, 2°). The error bars 
show the uncertainties of the individual MLS profiles. Approximate potential temperature is 
shown near the right vertical axis. Locations o| the matched MLS profiles are shown as well. 
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Figure 5: Vertical profiles of the difference between ER-2 and MLS ozone measurements in 
ppmv (panel a) and in % (panel b) using the THT with the match criterion of (2 h, 2 °, 2°). 
The error bars show only the random errors determined in Appendix. Numbers of matches are 
shown near the right vertical axis. 5 
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Figure 6: The same as in Figure 4, but for the POAM-III and ER-2 ozone measurements on 
February 2 and January 23, 2000. 
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Figure 9: The same as in Figure 7, but for the POAM-III and ER-2 H 2 0 measurements. Cyan 
(green) lines and symbols correspond to the Harvard (JPL) measurements. 
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Figure 10: The same as in Figure 4, but for the ER-2 and MLS CIO measurements. 
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Figure 11: The difference between the model-calculated and MLS-measured CIO for the forward 
trajectories from the ER-2 locations using the NAT (pink line) and STS (black line) PSC 
schemes. The AER box model calculations are constrained by the ER-2 measurements as 
mentioned in the text. 
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Figure 13: (panel a) The difference between the MLS and ER-2 gaseous HN0 3 measurements 
determined by all (red curve with errorbars), forward (blue curve), and backward (green curve) 
trajectories assuming that HN0 3 is a passive tracer, (panel b) The same as in panel (a), but 
for the forward trajectories only. The pink and black lines shows the model results with the 
NAT and STS PSC schemes, respectively. The yellow line shows model results with the STS 
scheme and zero intial content of HNO3 in aerosol. The blue line with errorbars is the same 
as in panel (a), (panel c) Profiles of the averaged temperature in the initial (labeled as ER-2) 
and final (labeled as MLS) points of the matched forward trajectories and the NAT threshold 
(labeled as NAT). 






